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Calculations of reaction rates for the third-order QED process of photon splitting 7 ~ ^ 77 
in strong magnetic fields traditionally have employed either the effective Lagrangian method or 
variants of Schwinger’s proper-time technique. Recently, Mentzel, Berg and Wunner [1] presented 
an alternative derivation via an S-matrix formulation in the Landau representation. Advantages of 
such a formulation include the ability to compute rates near pair resonances above pair threshold. 
This paper presents new developments of the Landau representation formalism as applied to photon 
splitting, providing significant advances beyond the work of [l] by summing over the spin quantum 
numbers of the electron propagators, and analytically integrating over the component of momentum 
of the intermediate states that is parallel to field. The ensuing tractable expressions for the scattering 
amplitudes are satisfyingly compact, and of an appearance familiar to S-matrix theory applications. 
Such developments can facilitate numerical computations of splitting considerably both below and 
above pair threshold. Specializations to two regimes of interest are obtained, namely the limit 
of highly supercritical fields and the domain where photon energies are far inferior to that for 
the threshold of single- photon pair creation. In particular, for the first time the low-frequency 
amplitudes are simply expressed in terms of the Gamma function, its integral and its derivatives. 
In addition, the equivalence of the asymptotic forms in these two domains to extant results from 
effective Lagrangian/proper-time formulations is demonstrated. 

12.20.Ds, 95.30. Cq, 97.60.Gb, 97.60.Jd, 98.70.Rz 


I. INTRODUCTION 


The third-order quantum electrodynamical process of photon splitting 7 — > 77 in a strong magnetic field, currently 
popular in several astrophysical models of different neutron star sources, was first studied over three decades ago. 
Due to analytic complexities encountered when investigating this interaction, it was not until the beginning of the 
1970s that a body of correct and uncontroversial results emerged. These early splitting calculations used either 
effective Lagrangian [2-4] or variations of Schwinger’s proper-time techniques [5-7], the expediency of which yielded 
compact analytic forms for the rates R when specializing to low energy (R oc w 5 ) or low field ( R oc B ) cases. 
After a hiatus of nearly two decades, photon splitting became of interest again in the literature [8-11] following the 
publication of an S-matrix calculation in the Landau representation of its rates by Mentzel, Berg and Wunner [1], 
specifically because of their contention that the earlier works cited above had seriously underestimated the strength 
of this process. The rates computed in [1] were later retracted in [12], with a sign error in their numerical coding 
having been discovered and corrected. Mentzel et- al.’s analytic derivation was the first comprehensive presentation 
of the application of a Landau representation technique specifically to magnetic photon splitting, though the QED 
formalism presented by Melrose and Parle [13,14] virtually provided an equivalent enunciation of such S-matrix forms 
for splitting amplitudes. More recently, W^eise, Baring & Melrose [11] confirmed the analytic derivation of [1]. The 
Landau representation calculations and most of the earlier effective Lagrangian and proper-time presentations were 
generally applicable to non-dispersive regimes below the pair creation threshold ( hu — 2mc 2 ), where the momentum 
vectors of the initial and final photons are collinear, and arbitrary field strengths. 

Below pair threshold, the effective Lagrangian approach of [2-4] and the proper-time calculations in [5-8] appear 
much more amenable for the purposes of numerical evaluation than the S-matrix formulation in the Landau repre- 
sentation. This arises because effective Lagrangian and proper-time (collectively referred to by the label ELP here) 
methods produce results that involve triple integrals over relatively simple (hyperbolic and exponential) functions, 
while the S-matrix amplitudes integrate over the parallel momentum p z and include a triple summation over the 
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Landau level quantum numbers of the intermediate pair states. Both techniques start from different but equivalent 
[15] forms of the electron propagator, and hence S-matrix computations [1,11] should yield identical results to proper- 
time numerics [2-8]. For the specific case of magnetic pair creation 7 -+ e + e” , such an equivalence of the S-matrix 
and proper-time methods has been demonstrated [16,17], but only via continuous asymptotic approximations that 
smoothly average out the exact “sawtooth” resonance structure. Vet the S-matrix Landau representation approach 
explicitly retains the resonances in the scattering amplitudes above pair threshold, whereas the ELP methods elimi- 
nate such information early during developments. Photon splitting becomes effectively first-order in a f at any one of 
a multitude of pair resonances, generated when the intermediate states become on-shell. Hence it is quite possible 
that splitting can compete effectively with pair creation as a photon absorption mechanism above pair threshold. As- 
certaining whether this is true is an interesting physics question. Moreover, if splitting is approximately as probable 
as pair creation above threshold, then it manifestly changes the character of vacuum dispersion, so that quadratic 
(and by inference perhaps higher order) contributions to the vacuum polarization tensor become significant relative 
to the standard linear ones used in the derivation [4] of kinematic selection rules for splitting. Hence, the generation 
of exact and compact expressions for the rates for 7 -> 77 valid both below and above pair creation threshold is 
clearly a worthwhile enterprise from a physics perspective. 

Developed expressions for the rates for photon splitting are also important for astrophysical applications of this 
process, particularly to effect efficient and accurate computations of such rates. These applications have so far focused 
on neutron star magnetospheres, primarily on models of soft gamma repeaters (SGRs) and strongly-magnetized 
pulsars, both being extremely topical in the astrophysics community at present. The potential importance of splitting 
in neutron star environments was suggested by [4,18,19]. Possible formation of splitting cascades has been explored 
in models of SGR transient outbursts as a means of softening the spectrum efficiently with no production of pairs 
[20-23]. If both polarizations can split, or if polarization switching is active during SGR outbursts, then the properties 
of the splitting cross-section guarantee emergent spectra in the observed range (20-150 keV) and of the observed shape 
for all fields in excess of around 10 14 Gauss [20-22], provided that the emission region is not concentrated near the 
polar cap. The spectral properties of SGRs in quiescent emission appear to be distinct from those during outburst. 
Pulsations and temporal increases of their periods (i.e. spin-down) have now been observed [24—26] for two of the four 
confirmed SGRs (SGR 1806-20 and SGR 1900+14), leading to inferences of fields in the vicinity of 10 Gauss. The 
connection between these pulsars of extremely high magnetization, so-called magnetars , and conventional radio/X- 
ray /gamma-ray pulsars is not well-understood. Baring & Harding [27] postulated that radio quiescence, a property of 
the SGRs, may be common in magnetars due to the efficient action of photon splitting and other effects in suppressing 
the creation of pairs. Photon splitting also has spectral implications for such pulsars with more modest fields. [28] 
demonstrated that the unusual absence of > 30 MeV emission in the gamma-ray pulsar PSR 1509-58 (whose spin- 
down field is ~ 3 x 10 13 Gauss) can naturally be explained by the operation of 7 -+ 77 in the intense magnetic and 
gravitational fields near its surface. 

Several desirable goals are immediately identifiable on the basis of this historical path for the study of the physics of 
photon splitting, and the needs of the astrophysics community. It would be satisfying (i) to obtain analytic expressions 
for rates that are valid above pair creation threshold using the Landau representation methodology, (ii) to know 
whether the analytic formalism of Mentzel et al. [l] can be developed and simplified, and (iii) to demonstrate a formal 
equivalence between this S-matrix Landau representation approach and extant results from proper-time/ effective 
Lagrangian techniques. This paper addresses these issues, using the verified analytic formalism of Mentzel et al. as 
the starting point for mathematical developments. The analysis here considers all the polarization modes that are 
permitted by the CP invariance symmetry ( J_->|||| , and ||— >-L|| ), and applies for collinear momenta of the 

incoming and outgoing photons, i.e. when the effects of vacuum dispersion are neglected. A significant development 
provided in this paper is the dramatic simplification incurred by algebraically performing the summation over the 
spin states that are incorporated in the electron propagators. The resulting expressions in Section II A (first stated in 
[11]) are relatively compact, and of an appearance familiar to Landau representation/S-matrix theory applications to 
magnetized environments (i.e. including associated Laguerre functions). Furthermore, here the integrations over the 
momentum parallel to the field are performed analytically for the first time in Section II B, rendering the splitting 
rates in most amenable forms (see Eqs. [12] and [15]) that are optimal for numerical applications: the analytic forms 
presented consist of just triple summations over Landau level quantum numbers of the intermediate states. These 
general results are valid both below and above pair threshold at non-resonant photon energies, and provide substantial 
advances over the work of [1]; they are much more suitable for numerical evaluation since many cancellations have 

been eliminated algebraically. . 

Two specializations are discussed in Section III, primarily to (partially) demonstrate equivalence of the Landau 
representation formalism presented here with extant proper- time/effective Lagrangian limiting forms for splitting 
rates, and simultaneously to serve as a check on the mathematical manipulations of this paper. Results are presented 
for all three polarization modes permitted by CP invariance in the limit of zero dispersion. The first asymptotic 
regime is (see Section III A) for highly supercritical fields, B » B c = m 2 c 3 /eh , where in the case of ±->|||| , the limit 
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was found to concur with a recent analytic result that was obtained by Baier et al. [ 8 ], while new results were obtained 
for the other two modes. In the second specialization, in Section III B, asymptotic results for energies u <$C me 2 well 
below pair creation threshold were obtained, reproducing the cubic energy dependence of the amplitudes obtained 
by other QED techniques. Moreover, new and compact expressions for the scattering amplitudes in this low energy 
limit are derived in terms of the logarithm of the Gamma function, its integral and their derivatives. These simplified 
forms in Eqs. (41) and (42) are also produced from extant integral forms for splitting matrix elements derived first in 
[2,4], thereby facilitating the first analytic demonstration of the equivalence of splitting rates obtained by the S-matrix 
formulation in the Landau representation and those derived using Schwinger- type techniques. 


II. THE GENERAL S-MATRIX FORMALISM 


The rates for photon splitting within an S-matrix formulation can be developed using a variety of conventions; 
here the Landau representation used by Mentzel, Berg and Wunner [1] is adopted, and formal developments lead to 
an independent confirmation of their analytic derivation. Specifically, for a field (0,0, B) , this approach uses a 
representation of the electron /positron wavefunctions sis eigenstates of the magnetic moment (or spin) operator \i 2 
(with = rrur -f 75 [3(T x [p + eA(x)]) in Cartesian coordinates within the confines of the Landau gauge A(x) = 
(0, Bx } 0) . Such states turn out to be very convenient because they generate useful symmetry properties; they were 
identified by Sokolov and Ternov [29], who dubbed them states of “transverse polarization. Let e , e and e denote 
the polarizations of the initial and final (primed) photons (e, e', e lt =J_, ||), and k ^ = (u;,k), k ^ = (uj ,k) and 
k n = (u/\k") denote the absorbed and produced photon four momenta. Then the total rate for splitting via the 

^ V 5 ' /g\ 

polarization mode e — > e f e f/ can be written, using Eqs. (27)-(29) of [1], in terms of the S-matrix element £y t - , which 
is the sum of six terms Sj^j corresponding to the six viable time-ordering possibilities: 


Re-*e'e n 


- f ^ ^ 2 mc2 

2 ^ 87T 3 ) h 


d 3 k' d 3 k" i 


E 

j — 1,6 


c(3) 


(i) 


where V and T denote the volume and time associated with the interaction calculation and the factor of 1/2 out 
the front avoids double counting of the final states. The priming convention adopted throughout the paper is one and 
two primes for the produced photons and no prime for the initial photon. Since the S-matrix element contains a delta 
function S A (k ^ — k ^ — fc") prescribing four-momentum conservation for splitting, it is squared in the usual way using 
|J 4 (it /i - k^ - fc")l 2 “> [V r T/( 27 r) 4 ] S 4 {k „ - k'p - *") . Note that the S-matrix element should possess a cubic dependence 
on photon energies when well below pair creation threshold, due to parity symmetry, photon gauge invariance, and 
the antisymmetric nature of the electromagnetic field tensor; details are discussed in [4]. 

Before writing down expressions for the S-matrix element terms, it is appropriate to identify the dimensionless 
convention that shall be adopted throughout this paper. Since the electron rest mass m is the only mass that enters 
into this QED problem, we opt to scale all energies by me 2 and momenta by me unless otherwise specified. This 
includes a scaling of mc 7 /k for photon frequencies u > . In the spirit of this convention, we choose to use the symbol 
e to represent, dimensionless electron energies and reserve E ( = emc 2 ) to denote “dimensional'’ energies as in [ 1 ]. 
In addition, the magnetic field will be expressed in terms of the quantum critical field B c — m 2 c 3 /eh hereafter, so 
that B — 1 denotes a field of 4.413 x 10 13 Gauss. 

The convention for polarizations is identical to that assumed in [1], who opted for real polarization vectors with zero 
time components. The polarization states JL and || are defined according to whether the photon’s electric vector lies 
either perpendicular or parallel (respectively) to the plane containing the photon’s momentum k and the (uniform) 
magnetic field B vectors, the convention of [1,13,6,16,10,30]. In the limit of zero dispersion, three polarization modes 
are permitted by charge/parity (CP) invariance in QED, namely ±->|||| , -L->±± and ||->_L|| . However, Adler [4] 
showed (see also [31]) that for weak vacuum dispersion (roughly delineated by B < 1 ), where the refractive indices 
for the polarization states are very close to unity, energy and momentum could simultaneously be conserved only 
for the splitting mode -L— H||| . This kinematic selection rule applies to gamma-ray pulsar magnetospheres where 
plasma dispersion is negligible. In magnetar models of soft gamma repeaters, where supercritical fields are employed, 
strong vacuum dispersion arises. In such a regime, it is not clear whether Adler’s selection rules still endure, since his 
linear dispersion analysis omits higher order (quadratic) contributions [13,14] to the vacuum polarization tensor (e.g. 
those that couple to photon absorption via splitting) that may become significant in supercritical fields. Furthermore, 
plasma dispersion effects, which can nullify the vacuum selection rules, may be quite pertinent [32] to soft gamma 
repeater magnetospheres, rendering them distinctly different from those of conventional pulsars. Therefore, in the 
interests of generality, consideration of all three CP-permitted splitting modes is adopted throughout this paper. 
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The derivation of the S-matrix element proceeds along lines identical to those in Mentzel, Berg k Wunner [1], with 
the result being an exact reproduction of their analytic formalism, as reported in Weise, Baring k Melrose [11]; for 
details, one is referred to [1], Of the six Sj^j contributions to Eq. (1), it is sufficient to explicitly present just one. 

* .5. .5. 4^ / *■ 

2>++(k W )^;»(k')^n*(k) - ^nn“(k")^ n (k')P n - + „(k) 


e e'e" (e + e" + w' — it) (s' + e" + w - it 1 ) 


( 2 ) 


I p'.=p.-K.p’.’=-P‘-*i 


where Of — e^/(hc) is the fine structure constant, and the energies £ and £o are defined in Eq. (4) below, with 
similar definitions for the primed energies involving primed quantum numbers and momenta of the virtual electrons. 
Using the J notation in Eq. (6) below, 

2>+,+ (k") = J{~K I n' - 1, | n, 0) K*k 4 + 4*k 2 ] e" 

+ •/(-*" | n', -k” | n - 1, 0) [4* /ci + k?k 3 ] e+ 

- | n', | n, 0) (4 *k 4 + 4'k 2 ] - | n' - 1, | n - 1, 0) [<k 3 + «i]}e? . 

V+-„(k') = | n - 1, 0 | n" , k' y ) [«K + 44"] 

+ J(-t; I n, 0 1 n" - 1, k' y ) [44" + 44 *] e' + (3) 

- { J(-*; I n, 0 I n", *> ) [44* + 441 - | n - 1, 0 I n" - 1, *;) [44" + 44"]}4 . 

X>~,t„'( k ) = J(k x | n" - 1, k' y | n', -fc") [<4 + 441 e - 

+ J(k x | n", 4 | n' - 1, -fc") [44 + 441 e+ 

- {J(k x | n", 4 | n', -4) [4«" + 44') - W* | n" - 1, 4 | n' - 1, -*") [44 + 44']}^ - 

where the polarization vector e ^ = (0,e x ,e y ,e z ) is specified by e± = t x ± ie y and e z , and similarly for the final 
photon polarizations (primed). The other three Ps in Eq. (2) are not displayed here for brevity; they can be obtained 
from those in Eq. (3) simply by the interchange e + *4 e_ of polarization, components (and similarly for primed 

components) and a relabelling of the J s that produces a correspondence U rl , n (k) — > n , (k) . 

Several notations need to be identified. First, the particles have energies £ , and momentum components p 2 along 
the field. The energies e and e 0 that appear here are, respectively, with and without the parallel momentum p z : 

2nB , (4) 


y/l + p? + 2 nB , £o = + 


with n denoting the Landau level quantum numbers, as usual. The other quantum number pertaining to the 
eigenstates of p z is a = ±1 , which signifies the spin state of the fermions ( [1] used the label t ; here the notation of 
Melrose and Parie [15] is preferred), and satisfies p z xp = <T£ 0 rp . It does not appear explicitly in e, but is embedded 
in the spinor coefficients : 


( Kl \ 






W 
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0 

1 


(5- s+ 
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\0 0 
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Pz y e + so 

\yi v/(e 0 - l)(e + eo) J 


( 5 ) 
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primed coefficients k' and /c" are similarly defined in terms of primed momenta and Landau level quantum numbers, 

subject to the momentum conservation implicit in Eq. (2). 

The J functions that appear in Eq. (3) are integrals over the oscillator functions (Hermite polynomial products), 
a form undeveloped in [1]. Here, Eq. (7.377) of [33] is employed to express these integrals analytically in terms of 
generalized Laguerre polynomials, L" " n (x) (see also Eq. (47) of [15]): 

J{cx\ri, 0) = exp(-z ^ [/? + /?']) j In '> n {w) ^ 


where p (> 0) and the phase ^ are introduced for convenience of notation: a = pcos^ and /?'-/? = psinip . 
Note that the fis are always k y s and a is always a k x . Here the / n ',n functions follow the Sokolov and Ternov 
convention [29] up to a factor of n \ , being related to the J functions of Melrose and Parle [15], and both are defined 
in terms of the generalized Laguerre polynomials (see [33]): 


7n\n(^) — (“I)™ ^n > n , { x ) — 



x (n' — n)/2 L n'-n (x ) 


n! > n . 


(?) 


Values for n > n* are obtained by interchanging indices, as indicated. Hereafter, the (modified) Sokolov and Ternov 
convention for writing the Laguerre polynomials will be adopted. Complex conjugation of Eq. (6) can be used to 
establish the identity 

SftOO = ( — l) n-n {^nn' _?, (k)}* , ( § ) 


noting that the k products are either purely imaginary or real, for all choices of spin quantum numbers. The three 
factors like (— l) n _n appearing in the second product of three Ds in Eq. (2) cancel, leading to this product being 
just the complex conjugate of the first three Ps, This useful symmetry property clearly underlines the convenience 
of the Sokolov and Ternov choice of wavefunctions when adopting real components for the photon polarization. 

The form of the contribution to the S-matrix element in Eq. (2) is identical to that for , given in Eq. (25) 

(3) 

of Mentzel, Berg and Wunner [1]. In the same fashion, it can be found that the expression derived here for S fi 2 
is absolutely identical to Eq. (26) of [1], thereby providing confirmation of their analytic developments, it can be 
obtained by using the substitutions k / ++ k /; , n f ++ n n and e f ++ e n in Eq. (3). All other contributions result 

from application of the cyclic permutations 


7 + i • e" -+ k'^, , k^ e* 

P b _v h* p' k f p' — y k” e" 

1 ■ ^ y 7 J 


v b” 

i 7 p 1 

k” p" -k 



( 9 ) 


where kp = (uq k) , = (0, e X} e y , e z ) , etc. Observe also that a minus sign and the complex conjugation of the po- 

larizations are always associated with the initial photon since it is absorbed in the process. Given these permutations, 
the crossing symmetry for splitting is manifested in the following relationship between the various terms like those in 
Eq. (2) that contribute to Eq. (1): 


qO) p q(3) 


q( 3) _ 

D Ji t 4 ~ 


p q( 3 ) 

^+1 ^/i, 1 


c(3) _ 


P + i s ( j% 


q( 3) _ 

t,6 


p c(3) 


(10) 


where the permutations act as operators. This symmetry can be expressed in a multitude of ways using the identities 

P+iP-i = /= F-iP+i and P|, = I. . 

It is important to remark that the derivation of analytic forms by Mentzel, Berg and Wunner is not the first in the 
literature relating to S-matrix applications to photon splitting. The papers by Melrose and Parle [15,13,14] dealing 
with various aspects of QED in strong magnetic fields, specifically from a wave dispersion/response tensor approach, 
constructed the S-matrix element for splitting in Eqs. (46) and (47) of [13], which incorporated the quadratic vacuum 
response tensor given in Eq. (36) of [14]. This tensor is obviously of a standard S-matrix Landau representation 
appearance. Eqs. (2) and (3) can be generated directly (and also ) from the Melrose and Parle evaluation after 
a modicum of algebra. Hence, Eqs. (2) and (3) here, and Eqs. (25) and (26) can be used as reliable starting points 
for further S-matrix developments. 
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A. Analytic Reduction: Summation over Spin States 


The form in Eqs. (2) and (3) is quite cumbersome. It can be simplified considerably by (i) specializing to specific 
but representative directions of photon propagation and (ii) analytically performing the summations over spin states 
<t , a' and a" . Restricting the photon motion to the x -direction yields photon motion perpendicular to the field: 
since splitting is collinear in the non-dispersive limit discussed earlier in this paper, it follows that. k z = k' z = k " = 0 . 
This choice dramatically simplifies coefficients of the Laguerre polynomials in Eq. (3). Without significant loss of 
generality, setting k y = k y = k” = 0 removes nearly all of the phase factors in the definition of the J s in Eq. (6), 

leaving just i n '~ n . Three such factors emerge in the triple product of V s, leading to a factor of (—1)" - 

The CP symmetry possessed by the splitting process becomes most evident at this point, since it is now simple to 
derive the CP selection rules. The specification of k y = k y = k” = 0 and k 2 = k' z = k" = 0 yields only one possible 
component of polarization perpendicular to the field, e± = £ y = — = and one conceivable component of 
polarization parallel to the field, e|| = £ z (and similarly for primed quantities). The polarization (electric field) 
vector of the photons is, of course, normal to the photon momentum vector, which automatically spawns the notation 
for the two possible polarization states: J_ : ej_ = 1, e|| = 0 and || : e± = 0, ey = 1 . Prom the presence 

of subtractions in the numerators of the integrands of Eq. (2) together with the complex conjugation property in 
Eq. (8) and the proportionality of the I>s to factors like i n '~ n , it follows that only terms with an odd number of 
ei factors contribute to , i.e. terms proportional to e.\\z\f'l , e \\ e 'i. e \\ > e - Le '\\ € \\ anc * ej - e i e i ■ All other terms 

cancel identically to zero. By virtue of the permutation symmetries in Eq. (9), this is also true for all other Sj t j . It 
is then trivial to deduce the OP selection rules for photon splitting, namely that the only permitted transitions are 

i~>ii , ±-H||| , ll-nil * ( ll ) 

The three other splitting transitions all have S-matrix elements that are exactly zero for collinear photon momenta, 
and hence are forbidden. This technique for CP selection rule derivation was implemented in [1]. These restrictions 
are simply consequences of the charge conjugation (C) and parity (P) symmetries of the splitting process, i.e. relating 
to the transformations k — > — k and B — > — B 

The summation over the spin states a, cr' and a" (= ±1 ) produces a dramatic simplification in the appearance 
of the S-matrix elements. Such spin summations act only on the products of the k, s that appear in Eq. (3); the 
algebra is lengthy but straightforward, being facilitated by pairing Sjfj terms with denominators that differ only in 
the sign of their photon energies. The total splitting rate in Eq. (1) can be written in the form 



where u" = w — u/ is implicitly understood from the conservation of four-momentum. While these rates will be 
expressed for photon propagation normal to the uniform magnetic field, the results for general photon obliquities 9 
to B can be obtained via a simple Lorentz transformation: u — > u> sin 6 , u >' —¥ ui' sin 9 , ui" — > ui" sin 0 , together with 
an extra multiplicative factor of sin# applied to the rate in Eq. (12). 

The momentum dependence in the integrands of Eq. (2) can be simplified by forming sums of the products of energy 
denominators. Separating such sums into real and imaginary parts via the representation , leads 

to the definition 
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for xj = 0 or 1 ; generalizations to complex (relevant to calculating splitting rates above pair threshold and near 

pair resonances) are routine. These manipulations yield the following compact forms for the M e -*e'e" coefficients in 


Eq. (12): 


M 


i-dlll - 


M 


i-tii = 


-j (-l)"''- n '|\/8n n'n"B 3 Z 0 A^ 1111 + \/2n"5 [j" + Zi - Zo] A 

n ,n f t n rr 

+VmB \j' - lx + Zo] Ag _+l111 + V2nB [j + h + Z 0 ] A^ 1111 1 

Y, (-l) n " _n '{V8 nnW'fl 3 Zo Af^ 1J - + \/2n"B [j" -1 1 -1 0 ] A 

n t n' ,n u ^ 

+VWB [j' + Zj +Z 0 ] Ag _+±1 + \/2^B [j + Z! +Z 0 ] A|" + - L - L | 


-l-IIII 

2 


2 J ” L 


^ |Ki|| = _! £ | V8nn'n"5 3 Z 0 A 1 ^ 1 " + VWB [j" + h - Z 0 


n,n } n 


A 


II-+J-II 


+%/2n'B J' + Zi + Zo] Af^l' + V^B J -1\ +Z 0 ] , 


(15) 


results that are to be used in conjunction with Eq. (12). The factor of —Bj 4 is introduced to render the scaled 
amplitudes positive, and also to afford a direct mapping onto limiting forms obtained [8,10] by the proper-time 
technique, as will become evident in Section III. The A^ cV ' are differences of triple products of generalized 
Laguerre polynomials (defined in Eq. [7]); for -L^|||| 


A 

A 

A 

A 


-L-HIII 

j 

-L-HIII 

2 


1- 

3 


_L- 

4 


bi — 1 ,n'— 1 

— I" 

— 'n-l.n'-l 
_ /" 

•*« — l,n # — 1 

— 1” 

-* n — 1 , n ' — 1 


I'n'^n In\n "- 1 — ^n,n' l,n-l Ai'-l.n" 

•fn' — l,n" ~ A»,n' ^n" t n — 1 

Ai,n' ^n",n 

Ifi n , n 4i'-l,n" “ ^n.n' ^n"-l,n-l j 


(16) 


where the Sokolov and Ternov representation of the associated Laguerre functions in Eq. (7) is used together with 
the priming notation 




V,n 


= /„ 



(17) 


thereby aiding brevity. For the J_— mode, 


^n",n- 1 Ai'.n"-! ” 




A 3 ^ 1X 


— \ ,n — l,n" Ai",n-1 ^n' p n"— 1 

= — ^n-l,n' ^n",n-l Ai'-l.n" 


A 4 — I nn t^ i j In' — \. t n n 1 y n' — 1 ,n — 1 j 


(18) 


and the results for the ||-»±|| mode are not explicitly stated since they can be obtained by exploiting crossing 
symmetries: the inverse of the permutation in Eq. (19) yields the transformation A x ^ — ► — A| ^ , A 2 ^ — > 
^IH-MI 5 -A!j' fi ' 11 , A^ 1111 A^ X|1 . The A-"* 6 '*” can alternatively be expressed using the J£ 

functions of Melrose and Parle as in [11]. Note that the potential subtlety of having to include factors of 1/2 for 
some contributions from ground intermediate states is eliminated by the specific choice of the Sokolov and Ternov 
wavefunctions. 

The comparative simplicity of the reduced form of the S-matrix element relative to Eq. (2) is both notable and 
comforting. Unlike Eqs. (25) and (26) of [1], this developed form of the splitting S-matrix element has an appearance 
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familiar to S-matrix applications of QED in the Landau representation to strongly-magnetized systems, with products 
of generalized Laguerre polynomials multiplied by simple combinations of energies and momentum components. Ex- 
amples of previous work bearing such familiar forms focus largely on lower-order QED processes and include studies 
of synchrotron radiation [29,30], single photon pair creation [17,29], and vacuum [34] and plasma [35] polarization. 

For the purposes of the analysis in the next section, it is pertinent to define the cyclic permutations 

ui — > —ui" , ui' —y —u> , u>" — > **>' , 
n -> n" , ri n , n" -+ n' , 


in the spirit of the P +1 permutation in Eq. (9). These permutations will appear repeatedly in the developments 
below, and lead to the following transformation properties of Eq. (13): 


E 


R 

- 1,-1 



- 1,1 





-E 


R 

- 1,-1 ' 


( 20 ) 


with E^ x being invariant, symmetries that are consequences of the arrangements of electron and positron propagators 
in the Feynman diagram for splitting. These translate into obvious mappings between 3 , 3 and 3 and an 
invariance of the /„ . It is also easily seen that under this cyclic permutation, the factor in braces in the summation 
for A4±_»n is invariant, while the equivalent factor in the summation for ,Ad||_>j_|| maps over (up to a minus sign) 
to the factor in braces in the A4j__t|||| summation. As will become evident in Section III, the remaining powers of 
— 1 in the summations do not provide any unsatisfactory interference in the limits of low photon energy (w <^C 1 ) 
and high fields (B > 1), so that permutation symmetry can be extended to the total amplitudes in these specific 
parameter regimes. 


B. Analytic Reduction: Integration over Parallel Momentum 


Further analytic development is not only possible, but also desirable, given that the integrations over the momentum 
p z parallel to the field can be expressed compactly in terms of elementary functions. Such tractability facilitates both 
numerical evaluations and the derivation of asymptotic limits. In proceeding, since results are sought at energies 
sufficiently remote from pair creation resonances, the imaginary parts of the denominators in the E A ,^ are dropped 
in all further considerations, i.e., we consider only the functions E^ = ReE Ai ^ . 

It turns out that carefully-constructed contour integrations in the complex p 2 plane do not facilitate the p z 
integrations. Hence the first step in integrating over p z is effected by the more cumbersome and less elegant approach 
of completing the squares and rationalizing the denominators using products of factors like (e' ± £ ± u>) . These 
factors define poles p l3 of the p- integration for i and j being some combination of n, n' and n" . Such poles fall 
into two types: pair creation ones (e.g. see [17]) that contribute only above pair threshold, due to the structure of the 
splitting rate, and cyclotronic ones that must be considered below pair threshold. The appearance of such cyclotromc 
poles is an artifact of the rationalization of denominators, so that they are really pseudo-poles of the subsequent 
analysis; a consistency check on the algebra is that the S-matrix element be effectively continuous across them. It is 
convenient to define energies that correspond to the p,j poles: 



AT = 1 + 2nB , Af' = 1 + 2 n'B , Af" = 1 + 2 n"B (22) 


is used for the purposes of abbreviation. Observe that, taking advantage of the subjectivity of such definitions, a 
minus sign appears in front of the expression for e„'„" , a choice that preserves symmetries induced by the mapping 
in Eq. (19) in the results that follow. These definitions spawn the following useful identities for the momentum poles: 


Pin' 

2 

Pn'n" 

2 

Pn n n 


= e 2 nn ,-M = e; 


2 

n*n 

2 


2 

nn 
2 

n'n 

2 n ~ = Z 


-Af' 

'n' — Af" 
,ii — Af , 


(23) 


which immediately imply the possibility of poles along the imaginary axis. In fact, p 2 nn> > -min{.V, Af'} , with 
equality for w" = \Af - A"| 1/2 , and likewise for the other poles. Note that for the one-vertex calculations of cyclotron 
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emission and single photon pair creation and annihilation, the requirement that such poles be real, corresponding to 
real components of particle momenta on external lines, is precisely what generates thresholds (e.g. [17]) and kinematic 
cutoffs (e.g. [30]) for transitions involving various states. 

The rationalization of the denominators yields relatively compact decompositions for these sums, after much can- 
cellation and simplification. They take the form 


yR 


= C\,u + 


2_ 

W 




e'e" + 1 




t"e) 


(24) 


the simplicity of which is contingent upon the energy-conservation restriction w" — w — a/ . Here 

W = wu /w" 4- uftf — wW - - 


(25) 


Identities such as W — — 2u;u/(£ n " n ' + £n"n) prove useful in the ensuing analysis. The c \ >#i and t\^ coefficients 
assume simple forms when expressed as partial fractions. Consider first the result for , which has the coefficients 


ci,i = 0 
c l , l 


jfr' _ £n"n , £ n"n' 

1 1 i = — 2 > 7 


Pn"n - Pz Pn'n" “ Pz 
7 


ie's" £nn” t €nn 

*i,i -15 72 + 72 


Pn u n Pz Pnn ' Pz 


*1,1 - ZJ 


^ + 72 72 


Pnn * Pz Pn'n 11 Pz 


(26) 


Observe that a cyclic symmetry is immediately apparent: is invariant under the permutation in Eq. (19), as is 

evident from its original definition in Eq. (13). Similarly, the algebraic developments yield coefficients for the 
sum, which appears in the ee f terms, as 


- St { gn i ,,n ' £n,n ' i : [e n »n' + *'] + € f n£nn ”i [£*"»-"] } 
1 ^ W \ pl. n „ -pV Pn"n Pz J 


*- 1,-1 


- 1,-1 


£n l, n l + U* 

‘72 — IV 

Pn’n" Pz 


(27) 


£nn' 

“5 7 

Pnn' Pz 


.e"e _ en"n-V 

— “7T5 — 


7 


Pn n n Pz 


&n'n 

“5 7 

Pnn ' - Pz 


The coefficients for the sum E R , , that appears in the ee"A 3 terms and the coefficients for the sum £ R _, that 
appears in the e'e"A 4 terms are similar: there is little need to state them explicitly, since the coefficients possess a 
relationship to each other due to the permutation symmetry enunciated in Eq. (20). 

Given these decompositions, it is now fairly straightforward to evaluate the integrations over p z , expressing them 
in terms of the an elementary function / with real arguments e,y : 


m , £) 



i 


sJP-N 

2 

v/V-£ 2 


; lo 6= 


Z + y/£ 2 - M 
Z - 7T 2 - V 


arctan 


{ } ’ 


if Z~ > Af , 
if 0 < Z 7 < M 


(28) 


for real € . The identity arctan z = ( 1 / 2 *) log e [(l + iz)/(l - iz)] with z = -Z/y/J7 - £ 2 has been used to map across 
the singularities at Z — i v A* (cyclotronic below pair threshold) and guarantee bounded and continuous behaviour 
of /(TV, Z)jZ at Z = 0 . The integral identity in Eq. (28) can be established quickly with the aid of result 3.513.2 in 
[33], using the substitution p 2 = VJ7 sinhf and partial fractions. Note that real values (either positive or negative) 
of Z are guaranteed by the formalism here, with Z = 0 being improbable due to the discreteness of the quantum 
numbers n , n' and n" . 

The integration of the coefficients I 0 of the A/ terms for each of the polarization modes are then straightforward, 
and the identities in Eq. (23) can be used to advantage. Similar terms appear in the T\ integrations of parts of the 
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coefficients of the other A; terms, which also possess integrands with terms proportional to 1/e , 1/e' and 1/e" that 
formally lead to divergences that cancel each other (an artifice introduced by the rationalization of the denominators). 
Using partial fractions, the divergent contributions can be written as integrals over the finite range -p < p z < p , 
rearranging to subtract off exactly-cancelling terms, and then taking the limit as p -> oo. Similar manipulations 
are used for the J" integration over where again the leading order terms are individually divergent yet 

collectively convergent. Partial fractions can again be used to enable rearrangements and separate the divergent 
terms, which are then integrated over finite ranges as with the X\ evaluation. The results are encapsulated in the 
identities 


1 o 

h 

J" 


— Jy -T i nn i + T, n'n" + ■F \ n"n j i 

= £ + yy '^Pnn' X- nn' + Pn'n" F r n'n" + Pn"n F n "n ^ i 

— £ — pp l Snn'Zn'n Fnn' + £n'n"{^n"n' + w n'n" + £nn"{^n"n 



(29) 


where 

nn' — /(A j ttin') T /(A' , £n'n)i 


Fn‘n" — /(A'', £n'n") + f(jV" , £n"n')t F n "n — /(A* > £n"n) + /(A , Snn"), 


(30) 


and 


£ = 


-77 log e N - 


-77 - log, A r ' - —T log. A/"' 

U) UJ UfUJ 


(31) 


No further integration is necessary: the cyclic permutations in Eq. (19) can be used to quickly derive expressions for 
J' and J from Eq. (29). 

At this point, it is salient to remark that the divergences at £ 2 = A r in the functions /( N,S) pose no problem 
for the integral evaluations in Eqs. (29), because these functions always appear two at a time. Below the pair 
threshold, these divergences are cyclotronic in nature, being encountered when ui — > |\/A ~ — vAV" | or for similar 
circumstances for the other photon energies. As u. tends to such a limit, for example, we observe that. 6 n'n" t \/A 
and £„»„» ->■ -sfW’ when M' > M" (without loss of generality). This opposition of signs guarantees cancellation of 
divergences when the arctan form of is used (arctan(l/z) — > 7t/2 — z as z — ¥ 0), so that continuity across 

cyclotron “pseudo-resonances” emerges naturally from Eq. (29), consistent with the continuity of the ^ functions. 
Continuity across pair resonances does not arise above pair threshold, so that true divergences emerge. 

The incorporation of Eq. (29) into the scaled matrix elements in Eq. (15) constitutes the final product of the general 
analytic developments in this paper, providing rates valid for all energies below pair threshold (and applicable for 
non-resonant energies above threshold), and for photon propagation normal to the uniform magnetic field. They are 
eminently suitable for numerical computations, having improved upon the analytic formalism of Mentzel, Berg and 
Wunner [1] (i.e. Eq. [2]) by performing the summations of the spin states and integration over the momenta parallel 
to the field that are associated with the electron propagators. Such developments are prudent prior to numerical 
evaluations due to the large degree of cancellation in these sums and integrations. 


III. ASYMPTOTIC LIMITS FOR HIGH B OR SMALL w 

A fruitful extension of this analysis is the exploration of the simplification of the scattering amplitudes and rates in 
two particular asymptotic regimes, namely the limit of highly supercritical fields, fl > 1 , and the specialization to 
photon energies well below threshold, i.e. w 1 . The benefits of such an investigation are twofold. First, it provides 
the first unequivocal analytic demonstration of the equivalence of splitting results from the S-matrix formulation in the 
Landau representation and effective Lagrangian/proper-time results from Schwinger-type formalisms in well-defined 
parameter regimes. In doing so, it serves as a powerful check on the developments here. Second, in the u C 1 case, 
it identifies a new, satisfyingly compact representation of the scattering amplitudes in terms of special functions that 

leads to an efficient means of computation. . 

These two parameter regimes are encompassed under the single limit lj 2 < 1+2 B , which thereby identifies 
the appropriate series expansion of the generalized Laguerre polynomials that appear in the amplitudes. For small 
arguments x , the leading order terms in the series for /„',„(x) can be found in the Appendix of [15]. Given that n , 
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n f and n n cluster in a manner such that \n f — n| ~ | n ,f — n\ ~ 1 , this series converges rapidly provided nx <§C 1 ■ 
Hence nw 2 /(2£) actually represents the true expansion parameter here, with u/ and a/' being similarly bounded. 
The leading order terms of such expansions for the Af -46 e are linear in the photon energies, while the next higher 
order terms are cubic; a more detailed exposition can be found in Weise, Baring and Melrose [11]. The series for the 
integrations of p z , namely 1 0 , 1\ , J " , J* and J (which do not depend on the polarization mode) are expansions 
in ur/( 1 -f 2 B) rather than w 2 /(2 B) . They are independent of photon energy to leading order, with a quadratic 
scaling with energy to next order. The series for lj 2 1 T 25 possess logarithmic character in the quantum numbers 
in situations when no two of them are equal (i.e. A 1 * ^ Af* ^ A f n ^ Af): 


77K > 


^ 4 \og p A r , 4 log fi AT 41og c A r 

Io ~ (A' - A") (A'" - /7) + (A" - h r ")(M -A 7 ) + (A'" - A') (A'' - A'") 

2AHok, AT 2A r 'log c A/ v 2A^" log. A'" 

Jl % K (N - N')(N" - A 1 ) W - M"){N - A') (A/"" - A") (A'' - A'") ’ 

and additionally involve inverse trigonometric functions when two ns (e.g. for N = A r ' ) are in fact equal: 


(32) 



II K — 


2X" 

{N - A r ") 2 


log e 


A' 2 
A” 77 A' - A"" 


+ 


4 A - - 

A'(A' - V") 



4 

(A* -AT) 



J" 


(33) 


where 


<?(*) 


arcsin j 
z\/l — x 2 


(34) 


and the identity arcsin x = arctan[x/\/l — x^] has been invoked. This retention of the inverse trigonometric functions 
is particularly relevant for determining the high B limiting forms of the scattering amplitudes. Relations similar to 
Eq. (33) exist for M = N l and Af f = Af" , obtained by the cyclic permutations through Afs and photon energies. 
The lengthier higher order (quadratic) terms are not explicitly stated for the sake of brevity. This concludes the 
preamble that guides the reader in the subsequent specializations. 


A. The Special Case of B > 1 


This regime is of particular relevance to the study of magnetars such as soft gamma repeaters. For the two modes 
j_^|||| anc [ ||-*_L|| , only the leading order terms for the A* and the momentum integrals presented in Eqs. (32) 
and (33) are required. Consider first the reduction of A4i~qm • Here the A 2 and A 3 terms contribute leading order 
terms only through n" = 1 , n - n f = 0 and ri = 1 , n == n" = 0 cases, respectively, where it is necessary to use the 
full forms in Equation (33), and inverse trigonometric functions appear through the Q(x) function, which assumes 
the arguments x = u ) f /2 and x — u/^/2 , A similar n = 1 , n r = n ,f — 0 term is identically equal to zero by virtue of 
the A 4 factor. The contributions from the Ai and A 4 terms possess an entirely different character, being infinite 
summations over n , with the values of n ; and n n being constrained by |n 7 — n| 4- \n u - n| < 1 , producing five 
groupings of the indices. The series is evaluated by truncating the sum at n < k , relabelling one of the logarithmic 
terms, and then taking the limit k — > 00 . The net result is (for u> < 2 ) 


M 





arcsin 


Of) 


+ 


4u;" 


u;V4-(u/') 2 


arcsin 




B » 1 


(35) 


which, when combined with Equation (12), yields the asymptotic high-B result derived by Baier et al. [8], and 
reproduced independently by Baring k Harding [10]; the overall rate for -L->-|||| approaches a value independent of 
B . Observe that the manifestations of the pair creation threshold for each of the final photons of || polarization 
(i.e. at u>' = 2 and u>" — 2) are the individually-divergent coefficients of the inverse trigonometric functions. Yet, 
collectively, due to the energy conservation relation w = ui' + u>" , such divergences cancel eac h other to yield a finite 
overall result as u — > 2 . For the incident photon of X polarization, the pair threshold of 1 + \/l + 2 B is remote from 
uj = 2 so that it would only become explicitly apparent when the amplitude was evaluated to higher order in B . 
Note also that the u <§: 1 limit of Eq. (35) is wu/w"/6 and reproduces results obtained in [4] and [6]. The functional 
form of Eq. (35) is plotted in Figure 1. 
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FIG. 1. The dependence of the scattering amplitudes for B » 1 , scaled by w 3 , on the fractional energy u//u : of one of 
the produced photons, for three different incident photon energies w (in units of me 2 ), as labelled. Only the two polarization 
modes with amplitudes asymptotically independent of B (in units of fl c ) in this ultra-quantum limit are depicted, namely (a) 
_L-*|| || and (b) ||->_L || ; their functional forms are given in Eqs. (35) and (36), respectively. The shape of the amplitude curves 
for 1-tlX is independent of w and is very close to that of the u = 0.1 curves in panels (a) and (b). While the ±-*|| || curves 
are necessarily symmetric about u / = w/2 , asymmetry is present in the ||-tl|| case where J represents the final photon of 
1 polarization. Note that the magnitude of diverges as pair threshold w 2 is approached. 


The equivalent result for the splitting mode M\\-*l\\ requires little additional algebra given that it can be obtained 
from the analysis just above using the cyclic symmetry transformations of Eq. (19). Carefully keeping track of signs 
and all photon frequencies by relabelling at the beginning of the manipulations, the roles of the A 4 and A 3 terms 
are interchanged, and the obvious result emerges: 


-Mu 


U) 


4u/' 


arcsin 


iO * 


infwl j — arcsin (-a-) + u/ 

V2 ) W 'V4-K) 2 \2) 


1 


(36) 


While not established before in the literature, the low energy limit of this, namely Af||->ij| « ww'w"/6 . yields the 
differential rate from previous expositions [4,6] of low energy approximations. The form of Eq. (36) is displayed in 
Figure 1, exhibiting the asymmetry expected under interchanges a/ u>" . In this case, pair threshold structure in the 
amplitude appears again for the two photons of parallel polarization (i.e. at w = 2 and oj" — 2), and is also absent 
for the produced 1 photon, being of higher order in B . Consequently, the amplitude possesses a real divergence at 
w - 2 , a noteworthy occurrence that is illustrated by comparing the two panels of Figure 1. Such divergences, which 
are not integrable over w (and therefore patently different in nature from the resonances encountered in rates for 
^ ), are characteristic of the photon splitting rate near resonances at and above the pair threshold of u i . 2 , 

corresponding to the creation of virtual pairs in various excited states. In fact, near such resonances, photon splitting 
necessarily becomes first order in o f like pair creation as the intermediate states go on-shell. 

The rapid increase of the rate of ||— >_L|| relative to that of T-t|||| is exhibited in Figure 2, where the rates have 
been scaled by the low energy (u> «; 1 ) limiting forms ( R{u) oc w 5 ) discussed in the next subsection. This particular 
scaling is chosen to illustrate deviations from the w « 1 asymptotic forms, and therefore to demonstrate the need 
for relinquishing use of them when sampling photon energies near pair threshold, a parameter regime very relevant to 
certain astrophysical calculations (e.g. see [27,28]). The dominance of the R\\->x\\ over R±^\\\\ near u -2 apparent 
in these B » 1 results becomes substantive in parameter regimes where the weakly-dispersive vacuum (i.e. for 
B < 1 ) polarization selection rules for splitting derived by Adler [4] (which prohibit ||->J-|| and 1-^-L-L splittings) 
may not apply if non-linear contributions to vacuum polarization or plasma effects are significant. This underlines 
the saliency of a detailed determination of the dispersive properties of the magnetized vacuum or plasma medium 
appropriate to a particular astrophysical scenario. . , , 

The derivation of the B > 1 form for the amplitude for _L->_L_L differs significantly from the results just expounded. 
First, contributions from = 1 , n = n' = 0 and n' = 1 , n = n" = 0 and n = 1 , n' = n" = 0 combinations are 
identically equal to zero by virtue of each of the associated A,- factors. This automatically implies that no inverse 
trigonometric functions that have arguments independent of B appear in the amplitude, a property not possessed by 
the other splitting modes. The consequences of this are twofold. First, this cancellation implies that the scattering for 
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±->_L_L is of a higher order in 5 than for the other two splitting modes. Second, since any potential appearance of 
inverse trigonometric functions spawned by the forms in Eq. (33) involves arguments that depend on B through the 
J\f s, these arguments are always small when B 1 , precipitating a redundancy with the low energy limit. Hence, it 
follows immediately that the scattering amplitude for J_— f-L-L in the regime of highly super-critical fields is identical 
to that of the B> 1 specialization of the low energy ( oj <§C 1 ) limit. As the latter has been derived in various papers 
in the literature (e.g. see [4,6,11] and the subsequent section), here it is sufficient to merely state the result. 


M 




WU Ul 

35 


5» 1 


(37) 


This extremely simple form differs profoundly from those of the other two modes because of the absence of photons 
of || polarization in the interaction. Hence any signatures of the pair threshold of 1 + \/l 4- 25 of _L photons are 
absent in the domain of cj < 2 , and a scaling-type form with obvious cyclic symmetry emerges. 



0 I i i^i — rT 1 1 — i — i — 1 T i I I I I 1 I — I — 1 

0 .5 1 1.5 2 


u 

FIG. 2. The total rates in the B » 1 limit for the modes 1-+|| || and ||-t-L||, computed according to Eq. (12) using 
the amplitude formulae in Eqs. (35) and (36), divided by the rates that would be computed when taking the low energy 
(u> C 1 ) limit of these amplitudes, i.e. A4x-»|||| ~ ww'w" /6 ss A4 j| _► x || • Deviations from such low energy approximations (i.e. 
R(w) oc u> 5 ), while significant for X-t|| || , are dramatic for ||-»± || near pair creation threshold w = 2 . 


B. Approximations for u> <C 1 

The low energy limit u> < 1 is of interest not only because it was the regime where compact analytic expressions 
for the splitting rates were first obtained [2-4], but also because the analysis that follows derives simple and ele- 
gant representations of the scattering amplitudes in terms of well-known special functions that provide a convenient 
alternative option for numerical evaluations. 

The amplitudes for each of the splitting modes should exhibit a cubic energy dependence [4] when u«l. Hence, 
a necessary product of the Landau representation formalism is that terms linear in photon energies should contribute 
exactly zero. For the polarization modes ±-f|||j and ||->--L|| , whose amplitudes are identical in the low energy limit 
[4,6], the demonstration of this is not dissimilar to the 5 » 1 analysis. The w«l restriction generates a single 
infinite series in n due to the clustering of n' and n" around n . The ensuing algebra in the simplification of this 
series is moderately lengthy, and requires re-indexing of the logarithmic terms to assume forms involving log e [l +2n5] , 
and also some relabelling of the rational functions. Care must be taken in these rearrangements due to the infinite 
nature of the series, and the technique adopted is outlined just below. The terms linear in photon energy result in 
zero, as expected: for more details, the reader is referred to [11]. The next order contribution is cubic in energy, 
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as desired, with terms coming from a mixture of (i) the linear terms of the A; combined with the quadratic higher 
order terms of the p, integrals, and (ii) the cubic A, terms in conjunction with the leading order (constant) terms 
from the p z integrals. The algebra is straightforward, but lengthy and tedious, generating an exact cancellation of 
all but terms proportional to ww'w" . This approach leads to a reproduction of the Cj listed in entirety in Appendix 
B of Weise, Baring & Melrose [11] for both the _L— >|||| and modes of splitting. Hence there is little point in 

replicating these expressions here; the reader is referred to [11] for details. 

These results are expressed as single infinite series in the label n , which sometimes starts at n — 0 , and sometimes 
begins at higher integer values (up to 3). Hence, an aesthetic goal is to rearrange some of these series so that the 
summations in each contribution begin at n = 0 , and then add the terms in the series together. This is a non- 
trivial exercise, given the divergent nature of the series in many of the individual contributions. Hence, considerable _ 
care must be taken when performing the rearrangements, for which there is no unique prescription. One choice for 
relabelling the sums is adopted by Weise, Baring & Melrose [11], though their end results expressed in their Appendix 
C do not facilitate analytic development in the most expedient manner, and were in fact, erroneous (discussed briefly 
below). An alternative and preferable choice for rearrangement, of the multitude of series over the label n is adopted 
here, outlined as follows. Inspection of the various Cj contributions in Appendix B of [11] reveals that they always 
consist of three types of terms: (i) logarithmic ones proportional to log e [l + 2(n + l)B ] , for l 0,il,i2,±3 , (ii) 
rational functions of 1+2 {n + l)B for t = 0, ±1, ±2 , and (iii) polynomials in n . A unique method for rearrangement 
is to truncate all series to finite ones with n < k , and then perform relabellings so that the first, two types of terms 
consist only of logjl + 2nB] terms and rational functions of 1 + 2 nB . This approach provides no particular focus on 
series that originate with labels n > 0 , but requires careful accounting of the remainder terms at the upper and lower 
ends of the sums, for which significant cancellation arises. The coefficients of the logarithmic functions, originally 
cubic in n , reduce to linear functions of n in this development. The consequent simplification of the series terms is 
counterbalanced by the transferral of complexity to the constant remainder terms, which are purely functions of k 
and B . Taking the limit of k -+ oo achieves the desired (and convergent) result. 

After considerable algebra collecting together all the constituent series in Appendix B of [11], and performing the 
rearrangement as just prescribed, one arrives at the following series representation of the scaled scattering amplitudes: 



for w -C 1 , where 

7T_*||||(n, B) 


Tj.-fii(n, B) 

defines the series terms. The remainders are quite lengthy, and are listed in Appendix A. 

Consider first the polarization mode X— 1|||| . While possibly only marginally simpler than Eq. (Cl) of [11], the 
series and remainder in Eqs. (38), (39) and (Al) naturally enable the development of a special function representation 
of the scattering amplitude. The finite summation over terms like ( x + n) log e (x + n) in Eq. (39) can be expressed 
using result 44.1.2 of [36] in terms of an integral of the logarithm of the Gamma function. At this juncture, the 
analysis begins to image parts of that generated in expressing the polarization properties of a magnetized vacuum 
via effective Lagrangian or proper-time techniques [37-39], as should be expected. Hence, it is appropriate to adopt 
definitions from such literature as much as possible. Following [37,38], here a definition for the generalized Gamma 
function Ti(x) of 

log e ri(z) = [ dt log e T(<) + i x(x - 1) - | log e 2 tt (40) 

Jo 

is adopted. Properties of this function, which include T^l) = 1 , are discussed at length in [40] and outlined in 

Rising Eq (Bl) one soon arrives at an expression for the scattering amplitude in terms of a handful of special 
(polygamma) functions, namely I^x) , and log e T(x) and its derivatives. This representation consists of two parts 
one independent of k , and one that involves a limit as k -+ oo of the remainder in Eq. (Al), combined with several 
terms incorporating the special functions with arguments that depend on k In evaluating this limit, most terms can 
be handled in a straightforward manner, and standard asymptotic series (e.g. see [33]) for log e r(x) and ip(x) as 
x -> oo prove useful. However, the treatment of the term involving the function log e Ti(l + k + 1/2 B) that appears 


+ logE [ n+ 2 £?] + (3 


1+2 nB 3(1 + 2nB)' J 


r 1 np n 4- 
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in the limit contribution is non-trivial, A series representation for this function for large arguments is required, and 
is presented in Eq. (B4). Assembling the various pieces, the limiting result as Ar -* oo is 



for u <C 1 • This is the sought-after compact analytic form that is comparable in simplicity to the one-loop effective 
Lagrangians calculated in [37,38]. Using series and asymptotic expansions for all the special functions present, it is 
routine to establish that M X -+m » (26B S /315) wwV' for B « 1 , while for B > 1 , one finds M X -> 1||| » , 

a result obtainable from Eq. (35). 

The developments are similar for the J_->±± mode: this representation again consists of two parts, one independent 
of k , and one that involves a limit as k — > oo of the remainder in Eq. (A2), combined with several terms incorporating 
polygamma functions with arguments that depend on k . This limit can easily be evaluated using asymptotic series 
to yield (for w 1 ) 


M±-> ii 


uuW> { ~W loge r (i) + 4F*(^) 

+ sF (A) + JS + W ~ W + if 7 ( loge 2?r + loge 25 ) } 


(42) 


Using series and asymptotic expansions for all the special functions present, it is routine to establish that .Mx-fX-L ~ 
(485 3 /315) uuj'w" for B < 1 , while for B » 1 , one finds « ww'o>'7(3 B ) , the result stated in Eq. (37). 

It must be remarked in passing that the expressions for J_— >-LJ_ in Eqs. (39), (A2) and (42) cannot be derived from 
the series in Eq. (C2) of [11], principally because that series expression is divergent, and therefore erroneous. Such an 
error was introduced by an inappropriate rearrangement of individually-divergent contributing series (leading to the 
addition of infinite contributions), a mistake that is avoided by the careful technique employed here in manipulating 
the results of Appendix B in [11]. Notwithstanding, the numerical results for the l->±± mode presented in [11] 
were effectively evaluated before any series rearrangement, and therefore remain valid. 



FIG. 3. The dependence of the scattering amplitudes, scaled by u>u'ui " , on the magnetic field (in units of 
B c = 4.413 X 10 13 Gauss) for the splitting modes J_-*|||| (see Eq. [41]) and 1-4 LL (see Eq. [(42]), for photon energies 
well below pair creation threshold. The amplitude for 1 1 — -L j| is identical to that for l - 1|| j| • 
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The compact analytic forms presented in Eqs. (41) and (42) represent the culmination of the w -C 1 focus here. It 
is the first time such simple forms for the scattering amplitudes involving just special functions have been calculated 
in this limit, though somewhat more convoluted, yet essentially equivalent, expressions have been put forward in 
[41]. A distinct advantage of the expressions in Eqs. (41) and (42) is the ease with which they can be accurately 
computed numerically. Their dependence on 0 is illustrated in Figure 3, replicating the numerics of [11] and earlier 
effective Lagrangian determinations [4], which are just as expedient since (see Eqs. [43] and [45] below) they involve 
just integrals of elementary functions. 

The low frequency result for .Mu^ill is not presented explicitly since it reproduces that for Adi-qui (e.g. see 
[4,6]); this is due to the crossing symmetries involved. Note that while the cubic dependences of all the modes at 
low energies reflect the lack of an energy scale in this domain (i.e. such as pair threshold), the normalizations are 
dependent on the polarization mode, particularly at highly supercritical fields where the M i-t-n amplitude is highly 
suppressed. This effectively represents how the rate normalization is sensitive to the (virtual) pair creation thresholds 
for the polarization states involved in a particular splitting mode. 

To conclude this presentation focusing on the w <C 1 specialization, an obvious objective is the re-derivation of 
Eqs. (41) and (42) starting with extant and well-known effective Lagrangian/proper-time (ELP) results, and thereby 
demonstrating analytically the equivalence of the S-matrix formulation in the Landau representation and Schwinger- 
type formalisms in the low energy limit. Consider first the mode _L — ^|j|j , for which such a determination is somewhat 
involved. The starting point is the integral expression [4,6,10] that corresponds to the scaled scattering amplitude 
that generates the same form for the rate as in Eq. (12): 


M 


ELP 

l-llll 


uu oj 


B 



f 3 s') cosh s 3 + 2s 2 s cosh s 

\ 4s + 6 ) sinh s + 12sinh 2 s + 2 sinh 3 s 


(43) 


which has and B » 1 limits matching those of Eq. (41). In the subsequent analysis, it is useful to manipulate 

integrations using the variable p = 1/B . The first step is to recognize that 1/s 3 times the factor in curly braces in 
Eq. (43) is a perfect derivative, namely dg/ds , where g(s) = { 1/ 4s 3 ) d[s coths — l]/ds — coth s/(6s) . Integration by 
parts is obviously the operative method, with the goal of retaining coths functions explicitly, combined with powers 
of s . After some algebra, one finds that 


mT p 




ds e 3 ^ b 


4 Ws 


_1_ 

B 



F 


1 _J _,J_ 

Tb? 60 + 3s 

r s 2i 

s coth s — 1 — g- 


(scoth 5—1) 


(44) 


results. The integral on the first line can be performed using identities 3.551.3 and 3.554.4 of [33], yielding the Gamma 
function and polygamma functions (or equivalently generalized Riemann Zeta functions) in addition to elementary 
functions. The only subtle part pertains to the second term of this integral, namely that contributed by the - 1/ (40s 2 ) 
factor. This can be differentiated with respect to 0 , evaluated to yield a ip function, and then the result integrated, 
noting the behaviour as 0 — ► 0 . The evaluation of the integral on the second line of Eq. (44) is much more involved. 
However, it has been performed before in the literature, and appears explicitly in calculations [37,38] of the one-loop 
effective Lagrangian describing refractive indices of the magnetized vacuum in QED. Hence the motivation for the 
particular partitioning of integrations chosen in Eq. (44). Details of the determination of this integral are found in 
Dittrich et al. [38], and the second line of Eq. (44) can be equated to —8 tt 2 /# 3 times the Lagrangian dS ( B ) (see 
Eqs. (2.4) and (3.16) of [38]), thereby introducing the Ti function. Collecting together the terms neatly generates 
an analytic form for that is identical to Eq. (41), so that the desired demonstration of equivalence of the 

Landau representation and effective Lagrangian forms is achieved. 

The procedure for the _L— ►_L_L mode is similar, though somewhat less involved. The equivalent scaled scattering 
amplitude obtained [4,6] from effective Lagrangian/proper-time techniques is 


M 


ELP 


uioj'ui" 

B 


f 00 ds —s/b f 3 cosh s ^ 3 — 4s" 
J 0 Y £ 1 4s sinh s + 4 sinh 2 s 


3s 2 | 

2 sinh 3 s J 


(45) 


Recognizing that the factor in curly braces can be written as -(3s/4) d[coths/s - l/s 2 )/ds + (s /4) d [coths - 
l/s]/ds 3 , integration by parts is again indicated, with identities 3.551.3 and 3.554.4 of [33] again proving useful. 
With manipulations similar to (but simpler than: the Fi function is not involved here) those for the -L-H||| , a 
modicum of algebra leads to the derivation of Eq. (42) from Eq. (45), as desired. This equivalence is a satisfying 
indication of the verity of the Landau representation analysis in this paper. 
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IV. CONCLUSION 


This paper has provided a detailed development of the S-matrix formulation of the QED process of magnetic 
photon splitting in the Landau representation, focusing on the case of zero dispersion where photon propagation is 
collinear. The formalism in Section II rederives and extends the exposition of Mentzel, Berg & Wunner [1]. The 
two principal general developments offered here are an analytic reduction via the summation over the spins of the 
intermediate pair states, discussed briefly in [11], and the analytic integration over the momenta parallel to the 
field incorporated in the electron propagators. This latter accomplishment is presented here for the first time. The 
cumulative product of these developments is a satisfyingly simple and elegant form in Eq. (15) for the scattering 
amplitude for each of the polarization modes permitted by CP invariance. These amplitudes possess products of 
generalized Laguerre polynomials that are common to QED processes in external magnetic fields, and elementary 
functions involving the photon energies and the various pair thresholds associated with the propagators. Moreover, 
the analytic forms presented consist of just triple summations over Landau level quantum numbers of the intermediate 
states, and are eminently suitable for accurate numerical computations both below and above pair creation threshold 
u> = 2 . The applicability of these results to regimes above pair threshold is a benefit of the S-matrix expansion in the 
Landau representation that is not afforded by effective Lagrangian and proper-time calculations: while these (latter) 
Schwinger-type techniques elegantly formulate splitting rates below pair threshold, they eliminate the resonance 
structure early on in their mathematical developments, a severe limitation above w = 2 . 

As an embellishment to these general results, specializations in two significant domains have been obtained. The 
first is for highly supercritical fields, B > 1 , reproducing in particular the result of [8] for the -L-H||| mode, and 
deriving new results for the other two modes permitted by CP invariance in the limit of zero dispersion. The second 
group of asymptotic results are for energies w < 1 well below pair creation threshold, where new and compact 
expressions for the scattering amplitudes have been derived in Eqs. (41) and (42) in terms of the logarithm of the 
Gamma function, its integral and their derivatives. These two domains of specialization herein have facilitated the 
first analytic demonstration of the equivalence of splitting rates obtained by the S-matrix formulation in the Landau 
representation and those derived using Schwinger-type effective Lagrangian/proper-time techniques. 
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APPENDIX A 


Here the remainders that appear in the series representation in Eq. (38) for the w < 1 specializations to the 
splitting amplitudes for polarization modes J_— >||[| and _L— are presented: 
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ix(*. B) = log e [* + 3 + 5^] 
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APPENDIX B 


In this Appendix, various useful properties of the r,(x) function, the integral of the logarithm of the Gamma 
function, that are needed in the u 1 specializations are stated. Given the definition of I\ in Eq. (40), it is 
elementary to establish , using 44.1.2 of [36], that 


k 

+ n) log e (x + n) = log e Fi(l + x + k) — log e Ti(x) . (Bl) 

n=0 

Taking successive derivatives with respect to x , one quickly arrives at well-known finite series representations of T(x) 
and its logarithmic derivative rf>(x ) ; see [33,42] for discussions of these functions and their series representations. 

An asymptotic series representation for the T\ function for large arguments is useful, and can be derived with the 
aid of the following series representation (see result 8.343.2 of [33]) for the logarithm of the Gamma function. 


1 1 °° Tfl 

log e r(x) = (x- |) log e *-*+ 5 log e 27T+2 (m+ lHm + 2) 


C(m + 1, x) - 


1 

~m+T 
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(B2) 


from which Stirling’s asymptotic expansion can be derived. Here, £(m,*) generalized Riemann Zeta function, 

defined in 9.511 and 9.521.1 of [33]. The integration of this series is effected using the identity C'( m >0 = -C(™+M) , 
and is mostly uneventful. However, the treatment of the m — 1 term in the summation is somewhat more subtle, due 
to the singular nature of C(0>0 > and requires taking the limit m — > 1 + , assuming m to be a continuous variable. 
Then result 8.362.1 of [33] comes in handy, and the series identity 


loge ri(ar) = ^ ( 2 * log e 1 1 - *) + J2 + i + 7E " l ] 


| ^ 

+ 2 51 (m + l)(m + 2) [' 


m= 2 


C(m, 1) — C("i, — 1 + 


(B3) 


follows, where 7 e = ~4’( 1) « 0.5772 is Euler’s constant. This series, which adequately substitutes for an asymptotic 
representation, can be used very effectively for numerical evaluations for all x > 1 . For the range 0 < x < 1 , this series 
also effects accurate evaluation of log, Ti(x) via use of the recurrence relation log, Ti(x) = log, 14(1 + x) - x log, x , 
an identity derivable from Eq. (40) with the aid of 6.441.3 in [33]. For large x , it then follows that 


log, r,(x) 


x(x — 1) 

2 


+ n) l0 Se^-T + il+ °(i) 


(B4) 


18 


where 


L i 



= lim 

k—too 
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Y1 n lo se n - 
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log e k + 



(B5) 


with numerical value L\ ~ 0.24875 . This is just the constant appearing in the magnetized vacuum polarization 
analyses of [37,38], where the Raabe integral form for it can be found. The second definition of L\ in Eq. (B5) can 
be obtained by setting x = 0 in Eq. (Bl), and is a result noted by [41], 
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